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The magnetic properties and critical behavior of both ferromagnetic pure and metallic nanopar- 
ticles having concurrently atomic disorder, dilution and competing interactions, are studied in the 
framework of an Ising model. We have used both the free energy variational principle based on 
the Bogoliubov inequality and Monte Carlo simulation. As a case of study for random diluted 
nanoparticles we have considered the Feo.5Mno.1Alo. 4 alloy characterized for exhibiting, under bulk 
conditions, low temperature reentrant spin glass (RSG) behavior and for which experimental and 
simulation results are available. Our results allow concluding that the variational model is suc- 
cessful in reproducing features of the particle size dependence of the Curie temperature for both 
pure and random diluted particles. In this last case, low temperature magnetization reduction was 
consistent with the same type of RSG behavior observed in bulk in accordance with the Almeida- 
Thouless line at low fields. A linear dependence of the freezing temperature with the reciprocal of 
the particle diameter was also obtained. Computation of the correlation length critical exponent for 
random diluted nanoparticles yielded the values z/=0.926±0.004 via Bogoliubov and i/=0.71±0.04 
via Monte Carlo. Differences are attributed to the so-called pairs-approximation in the variational 
model. From both approaches, differences in the v exponent of Feo.5Mno.1Alo. 4 nanoparticles with 
respect to that of the pure Ising model agree with Harris and Fisher arguments. 

PACS numbers: 75.10.-b, 75.50.Lk, 75.40.Cx, 75.40.Mg, 75.50.Bb,75.50.Tt 



I. INTRODUCTION 

The magnetic properties of bulk metallic systems hav- 
ing concurrently atomic disorder, dilution, competing 
interactions and characterized for exhibiting spin glass 
(SG) behavior have been widely studied from different 
points of view: experiment, theory and numerical simula- 
tion [lJ-Q . Among the systems with such characteristics, 
we can mention for instance CuMn B, F eAu @, FeAl 
[141], FeNiMn [13, [HI and FeMnAl [IlRi|] alloys. These 
alloys are interesting due to the richness of magnetic 
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phases that can be found such as ferromagnetic, anti- 
ferromagnetic, superparamagnetic, cluster glass, SG and 
RSG depending on stoichiometry, microstructure, degree 
of dilution, atomic disorder, magnetic field and tempera- 
ture. Typical I sing spin glass systems like those based on 
FeMnTiO, [lfllig, Fe(Cu i Al)Dy JS, LiHoYF [ljj, FeAl 
ESSE! and FeMnAl [H|l2>[l3 are S ood candidates 
to study SG and RSG related properties through Ising- 
based theoretical models where good agreement with ex- 
perimental results has been achieved. In particular, pure 
SG and RSG behaviors, in ternary FeMnAl alloys, arise 
from several ingredients including random atomic distri- 
bution of the alloy constituent elements in the crystalline 
structure, dilution provided by Al atoms giving rise to 
bond randomness and, finally, competition among the 
different exchange integrals involved. On this last re- 
spect, competition is given, essentially, by the difference 
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in sign and magnitude of the Jwe-Fe, Jpe-Mn and JMn-Mn 
exchange integrals. For high enough iron contents, a RSG 
behavior within the ferromagnetic phase governed by the 
Fc matrix can arise (l3j . 

Up to date, works reported on these kind of alloys deal 
with magnetic properties under bulk conditions. How- 
ever, to our knowledge, and despite of all the literature 
related with the so-called surface spin glass-like behav- 
ior in nanoparticles, no studies on metallic nanoparti- 
cles having concurrently bond competition, magnetic di- 
lution and atomic disorder within the entire volume of 
the nanoparticles have been reported. This fact has led 
us to consider the interplay between these effects and 
those arising from finite size when considering nanopar- 
ticles having such ingredients. The purpose of this article 
is to characterize from the magnetic standpoint, how the 
SG behaviors, found in systems like FeMnAl alloys under 
bulk conditions, become revealed in nanoparticles where 
the surface to volume ratio becomes increasingly impor- 
tant. Both the free energy variational method based on 
the Bogoliubov inequality and a Metropolis Monte Carlo 
simulation in the framework of a nearest neighbor Ising 
model were considered. The former, as energy minimiza- 
tion tool, has been already successfully employed in de- 
scribing the magnetic properties of this kind of systems 
where theoretical magnetic phase diagrams are in good 
agreement with the experimental ones [III EH EH • The 
layout of the paper is as follows. In Sec. II we describe 
the theoretical model and we emphasize the importance 
of a relationship for the average nearest neighbor coor- 
dination number as function of the particle size. In Sec. 
Ill we present our numerical results. This section pro- 
vides finite-size scaling analysis of both pure ferromag- 
netic nanoparticles, with an application to Ni nanostruc- 
tures, and Feo.5Mno.1Alo. 4 nanoparticles from both ap- 
proaches. Conclusions are finally presented in Sec. IV. 

II. THEORETICAL MODEL 

Several features lead us to consider an Ising model: i) 
magnetic frustration can be better resolved with an Ising 
model than, for instance, with continuous spin models [![, 
ii) it is in agreement with the framework we are interested 
in, which consists of iron-based nanoparticles with a very 
high effective magnetocrystalline anisotropy and where, 
despite of having cubic structure, a single easy axis can 
be experimentally induced |2l| - t23| . iii) it has been al- 
ready used in similar systems and quite good agreement 
with experimental data (magnctomctric measurements 
and hyperfine fields from Mossbauer spectroscopy) has 
been achieved 0, El El El, and finally iv) it allows to 
keep computational requirements under reasonable lim- 
its. Thus, our model is based on the following TV spins 
Ising Hamiltonian: 

N 

H = - Jjj(Ti<Jj -hy /Tj. (1) 



The first sum runs over nearest neighbors and <ii 

takes on the values ±1 or depending on whether the 
i th site is occupied by a magnetic atom (Fe,Mn) or an 
aluminum one, respectively. The exchange integral Jjj 
obeys the following probability distribution function ac- 
counting for disorder and the different couplings involved 

P(Jij) = p 2 S(Jij - Jfc-Fc) + 2pxS(Jij - Jfo-Mii) (2) 

+ X 2 6{Jij - JMn-Mn) + (<? 2 + 2xq + 2pq)S(Jij), 

where p, x and q, with p + x + q = 1, are the frac- 
tional concentrations of Fe, Mn and Al atoms, respec- 
tively. The terms p 2 , 2px and x 2 are the probabilities 
of having nearest neighbors Fe-Fe, Fe-Mn and Mn-Mn 
bonds, respectively, and interacting through the corre- 
sponding exchange integrals Jf c -Fc, Jfc-Mh and JMn-Mn- 
Here, Jfc-Fc hereafter simply J, was set to 12.846 meV 
only for pure iron nanoparticles with body-centered cu- 
bic (bec) structure, and was set to 16.872 meV for 
Feo.5Mno.1Alo. 4 nanoparticles having the same structure. 
These values reproduce the Curie critical temperatures 
of the corresponding systems under bulk conditions and 
the difference among them is attributed to the presence 
of both Mn and Al atoms and to the fact that Fe and 
Feo.5Mno.1Alo. 4 have different lattice parameters [l4| . 
Additionally, for the alloy, the remaining exchange in- 
tegrals were set to JFe-Mn=— otJ and JMn-Mn=— AJ with 
a=0.005 and A =0.03, [HJ and they correspond to the so- 
called competitive parameters. In this work, such values 
are kept fixed regardless the size of the nanoparticles to 
be considered. The last coefficient q 2 + 2xq + 2pq stands 
for diluted bonds, with Jy=0, corresponding to nearest 
neighbors Al-Al, Al-Mn and Al-Fe pairs. Finally, the 
second term in Eq.([T]) is the Zeeman contribution deal- 
ing with the coupling of the spins with a uniform external 
applied magnetic field h. 

Following the ideas proposed by Ferreira et al. [24| , in the 
pairs-approximation, the system is considered as formed 
by n\ single spins (S) and ri2 linked pairs (P) of spins with 
a total number of spins N = n\ + 2n2- Additionally, it is 
assumed that the magnetization can be obtained either 
from single spins or spins belonging to a pair. Thus, the 
trial Hamiltonian can be written as: 

Ua = -Is X! \- J ik°jVk + lp (crj + (Tfe)] , (3) 

ieS i.feeP 

where 7 S and j p are variational parameters which can 
be interpreted as molecular fields to be determined from 
energy minimization conditions. Here, the first sum 
runs over single spins and the second sum runs over 
pairs. Both Hamiltonians, Eqs.(fT]) and (J3|), can be re- 
lated through the variational approach based on the Bo- 
goliubov inequality: 

[F] < [F ] + [(H - H ) ] = m , (4) 

where F is the Helmholtz free energy defined by H, Fq 
is the free energy defined by Ho; (...)o refers to the ther- 
mal average in the ensemble defined by 'Ho, whereas [...] 
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represents a configurational average in which atomic dis- 
order is considered. According to the way as the system 
has been figured out we have: 



F = -k B TlnZ Q = -k B Tln(Z?- 2n *Z;z), (5) 



where 



2cos/i(/?7 s ), 



Z p = 2e 0J ^cosh(2p lp ) + 2e" /3J « 



(6) 



(7) 



are the trial partition functions for single and pair spins, 
respectively. The configurational average of any observ- 
able A is obtained from: 



[A}= I 



AP{Jij)dJ, 



(8) 



Calculation of the quantities [F ], [{H — ~Ho)q] and [$], 
following the same procedure as it has been described 
elsewhere [H, [H[ j leads to the following expression for 
energy minimization: 



dm 



2 {n 1 — n 2 ) (p 2 — 2pxa — Xx 2 )Jm 

Nh + (N - 2n 2 ) ls + 2n 2lp = 0, (9) 



where n' is the number of nearest neighbors, which de- 
pends on the crystalline structure, the type of boundary 
conditions and the system size. As concerns to nanopar- 
ticles we consider free boundary conditions. Since [$] 
diminishes as n 2 increases, we take n 2 as large as physi- 
cally possible, i.e. n' — n 2 . Thus, the number of linked 
pairs is maximized. Hence we obtain the following rela- 
tionship between the variational parameters or molecular 
fields 7 S and j p : 



Is 



{2 lp n'/N) - h 
(2n'/N) - 1 



(10) 



Our system is a spherical nanoparticle composed by N 
atoms arranged in a bcc structure with core coordination 
number z=8 and interatomic spacing a, which in the case 
of pure iron is around 2.86A whereas for Feo.5Mno.1Alo. 4 
is around 2.96A[14j]. On the basis of such representation 
and in order to get an expression for the maximum num- 
ber of pairs divided by the total number of atoms, i.e. 
n' /N, or, analogously an effective coordination number 
z c ff = 2n'/N, we have simulated particles with bcc struc- 
ture and different diameters D, in units of the lattice 
parameter a, and have counted the number of nearest 
neighbors pairs. Figure 1 shows the size dependence of 
z e g from which the following relationship is fulfilled for 
D > 3.4: 



ZcS 



2n' 



b 

D' 



(11) 



where z=8 and the best fit yields b — 10.13 ± 0.04 for a 
bcc lattice. In the case of a face-centered cubic (fee) 



6 - 



"S3 4 



bulk (bcc) 



O Eq.(17) 
Eq.(ll) 




0.0 0.1 0.2 0.3 

l/ P . 



10 12 14 16 18 20 

D 



FIG. 1: Size dependence of the effective (Eq. (|ll[) ') and the av- 
erage (Eq. (|17p ) coordination number for nanoparticles having 
bcc structure. Particle diameter D is given in units of the lat- 
tice parameter a. 



lattice, the following values must be used, z=12 and 
b = 12.65 ± 0.06. In principle, z c g can be interpreted 
as an effective coordination number for nanoparticles of 
diameter D whose core coordination number is z. Thus, 
the relationship between the molecular fields 7 S and 7 P 
can be rewritten as: 



7s 



IpZcS - 



z c s - 1 



(12) 



Magnetization can be computed either from single spins 
or from spins linked to a pair, and it must be the same: 



1 dlnZ, 



1 dlnZn 



(3 d ls 2/3 8 lp 



(13) 



After calculating the derivatives and performing the con- 
figurational averages using the bonds probability distri- 
bution function given by Eq. @ , we obtain the following 
transcendental equation for the magnetization: 



tanh((3j s ) 
sinh(2f3~f p ) 



cosh(2(3 lp ) + e- 2 P J 



2px 



cosh(20j p ) + e 2a P J 



2pq + 2qx 



cosh{2j3 lp ) + e 2X ? J cosh{2j3 lp ) + 1 



(14) 



Roots of this equation were obtained by using the Find- 
Root tool of Mathematica™. It must be stressed that 
such magnetization, according to the bonds distribution 
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function in Eq. ([2]) , corresponds to an average magnetiza- 
tion per bond whereas for the pure case (p=l) becomes 
a magnetization per site. The presence of crossed terms 
involving the atomic concentrations of the constituent 
elements p, x and q in Eq.(fT4f reflects the average over 
all possible nearest neighbors pairs as well as the random 
atomic distribution feature. Zero field magnetic suscepti- 
bility was obtained according to \ — (dm/dh) Q , yielding: 



X = 



^) ={(l-z cS )(p S ech*( ( 3 ls )y 1 



+ ^ [tanh{P ls )coth{2p lp ) 



sinh 2 {2jj-f p ) 



(cosh{2(3 lp ) + e- 2 /") 2 



2px 



{cosh(2j3 lp ) + e 2a ? J ) 2 {cosh{2j3 lp ) + e 2X P J ) 2 



+ 



q 2 + 2pq + 2qx 
(cosh{2f3 lp ) + If 



(15) 



The calculation of the Curie temperature Tc from Eq. 
(|14[) . for which m=0, is performed by taking the limits 
7 p —5-0 and j s — >0. This yields the following expression 
for the magnetic phase diagram: 



ZeS 



2px 



2(z c ff - 1) 



1 + e-WcJ 
x 2 

1 _|_ e 2A/3 c J 



I _|_ e 2a/3 c J 



2pq + 2qx 



(16) 



Here, /3 C = (fcsTc)- 1 . 



III. RESULTS AND DISCUSSION 



A. Coordination number 



where rii is the number of atoms having coordination 
number (see Fig. 2). Hence, by comparing the results 
derived from Eqs. (11) and (17) in Fig.l, we conclude 
that what we have called an effective coordination num- 
ber can indeed be considered as an average coordination 
number, i.e. z e s — (z). For diameters below around 
3.4a the average coordination number is characterized 
by jumps for which the discrete character of the system 
becomes more evident. A typical particle with D=1Q is 
illustrated in Fig. 2, where surface atoms with differ- 
ent coordinations are depicted with different colors. Our 
results are in agreement with those reported in metal 
Pt nanoparticles with different diameters where coordi- 
nation numbers for the first through fifth coordination 
shells were obtained by extended X-ray absorption fine 
structure (EXAFS) spectroscopy (25| . 





r 



FIG. 2: (Color online) Particle with diameter 75=10 having 
1067 atoms, bcc structure, and an average coordination num- 
ber (z)=6.96. Surface atoms with 4, 5, 6 and 7 nearest neigh- 
bors are colored green, blue, red, and yellow, respectively. 



From the particle size dependence of z e ff plotted in 
Fig. 1, particles with diameter D=10 (around 3 nm) al- 
ready exhibit an effective coordination number around 7, 
which corresponds roughly to 88% of that of the bulk. 
Above 13=10, the effective coordination number resem- 
bles that of the system under bulk conditions. Below 
that value, coordination number decreases rapidly and, 
therefore, strong modifications on the magnetic proper- 
ties are expected to occur in this range, i.e. below around 
3 nm. In order to gain a deeper insight on the interpreta- 
tion of z e g, we have computed the coordination number 
per particle Zi by counting the number of nearest neigh- 
bors surrounding the atom at the i th position, i.e. within 
the first coordination shell, and an average coordination 
number was computed according to: 



1 N 

8=1 



(17) 



B. Pure nanoparticles 

The particular case of pure ferromagnetic nanoparti- 
cles is easily obtained by setting p=l, q=0 and x=0 in 



k B T c (D) 



2J 



In [z e ff/ (z cS - 2)] ' 



(18) 



which gives the particle size dependence of the Curie 
temperature for nanoparticles with core coordination z, 
diameter D and a nearest neighbors exchange integral 
J. We want to stress that these results are not exclu- 
sively applicable to pure cv-Fe nanoparticles and they 
can be, in principle, employed for other pure ferromag- 
netic nanoparticles like Ni or Co. We can also relate 
this critical temperature with that of the system under 
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bulk conditions [lz 
temperature: 



in order to obtain a reduced critical 



T C {D) _ J(D)ln[z/(z-2)] 



T c (oo) J(oo)ln[z e g/(z eS -2)Y 



(19) 



where we have assumed that the exchange integral (J — 
J(D)) in the nanoparticle can be different from that un- 
der bulk conditions (J(oo)). In a first order of approx- 
imation, and by assuming the same exchange integral 
value, which could be reasonable for high enough parti- 
cle sizes, we have: 



Tc(D) 
T c (oo) 



Zeff 

z 



(20) 



if we assume that no structural transition occurs as a 
consequence of reducing the size. Otherwise, different 
core coordination numbers should be considered and the 
model is still applicable. Figure 3 shows the reduced 
critical temperature for different diameters according to 
our model. A comparison between Eq. (fT9|) for bcc (z=8) 
and fee (z=12) lattice structures and the approximate 
expression given by Eq.(f2"Uf is included. 
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FIG. 3: Particle size dependence of the reduced critical tem- 
perature for pure ferromagnetic nanoparticles with 2=8 and 
z=12, according to Eqs. (fT9|) and (f20)l . 

As is observed, Curie temperature (Tc) decreases as 
the particle becomes smaller due basically to the decrease 
in the magnetic bonds density. Therefore, the energy cost 
to carry out the transition is lower, and thus the critical 
temperature is also smaller. As D increases the critical 
temperature Tc tends to its bulk value. Other models 
have been already proposed to understand the mecha- 
nisms lying on the effect of the breaking of exchange 
bonds upon the Tc(D) function for nanoparticles. On 



this respect, the following expression has been proposed 
321 : 



Tc(D) 
Tc (00) 



3AL 
2D ' 



(21) 



where AL is the thickness of surface layer, and it has 
been considered as a parameter to characterize the de- 
ficiency in the number of exchange bonds for atoms at 
the surface region of a nanoparticle. However, this model 
can not reproduce successfully experimental data of mag- 
netite nanoparticles of different size with a constant AL, 
and hence, it has been suggested that AL should vary 
with the particle size, but such dependence has not been 
yet established 32] . Another model, based on the energy- 
equilibrium criterion between the spin-spin exchange in- 
teractions and the thermal vibration energy of atoms at 
the transition temperature and a size-dependent Debye 
temperature function, was developed in order to obtain 
both Tc{D) and Tj^(D) of ferromagnetic ad antifcrro- 
magnetic nanocrystals. Such model yielded the following 
expression for nanoparticles (27j : 



Tc(D) 
T c (oo) 



= exp[-(a-l)/(D/D -l)}, (22) 



where a is a measure of the root-mean-square (rms) ther- 
mal average amplitude of surface atoms vibration relative 
to the core and Dq denotes a critical size at which all 
atoms of the nanocrystal are located on its surface. Dif- 
ferently, our model contains just one adjustable parame- 
ter (J). Concerning a comparison with experimental re- 
sults, we want to stress that, in general, is rather difficult 
to obtain a diameter dependence of the Curie tempera- 
ture due to several factors like shape inhomogeneities, 
size distribution, and, in some cases like Fe nanoparti- 
cles, surface oxidation (26|. Despite of that, our results 
are, qualitatively, in good a gree ment with some others 
reported for nanostructures [27H30J. In order to evalu- 
ate the reliability of our model, we have carried out a 
comparison with some experimental data available for 
Ni nanostructures 0, |3l|. To do this, we have employed 
Eq. (IT5)) with z=12 corresponding to a fee lattice accord- 
ing to the crystalline structure of Ni, Tc(co)=631 K, a 
lattice parameter et=3.52 A, and we have also assumed 
that J=J(oo). Results are shown in Fig. 4 where we have 
also included the results from the models cited above. 

As observed, the agreement is rather good despite 
of the simplicity of our model and without considering 
any adjustable free parameter. Discrepancies can be at- 
tributed to the fact that, firstly, we have considered a 
simple nearest neighbors Ising model. Second, experi- 
mental data correspond to Ni nanostructures that are 
not spherical at all, whereas our model has been devel- 
oped for spherical nanoparticles. Third, for real nanos- 
tructures, the average lattice parameter is certainly dif- 
ferent from bulk, mainly for very small particles of some 
few nanometers of diameter and, therefore, exchange in- 
tegral should undergo changes. Such changes can also 
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shows a log- log plot of the reduced temperature (Tc(oo) — 
Tc(D)) /Tc(oo) versus particle diameter D, illustrating 
that the data obtained from Eq. ([TO)) follow the finite-size 
scaling relation (4p| - |43 |: 



T C (^)-T C (D) 
T c (oo) 



aD 



-x/v 



(24) 



from which our best estimate for the critical exponent as- 
sociated to the correlation length is v = 1.0001 ± 0.0001 
for both bcc and fee lattices. The observed exponent 
is slightly lower than the re por ted experimentally for 
Ni nanostructures (^=1.064)[31j, and very similar to 
that of a two-dimensional (2D) Ising model (f=l) but 
much greater than the observed in the three-dimensional 
(3D) Ising model (i/=0.6289)[Hj and mean field theory 
(^=0.5), suggesting a 3D— i-2D dimensionality crossover. 



FIG. 4: Semi-log plots of the dependence of the Curie tem- 
perature on diameter for Ni nanostructures. Comparison be- 
tween our model predictions (black solid line) without ad- 
justable parameters according to Eq. (|19[) . available experi- 
mental results for Ni nanoi ^articles (circles, Refs.[27|, HJ), Ni 
nanorods (squares, Refs.[27], (HI) and the theoretical mod- 
els described in the text with AL=0.8084 nm (dashed line, 
Eq.flHJ, Refs. H3,j3l) and D =1.4952 nm, a=1.811 (dotted 
line, Eq.((22}, Ref.[27[). Gray line stands for our model using 
Eq.JH) with a=0.6. 



be induced by other facts like volume magnetostriction 
of Ni. Regarding a particle size dependence of an effec- 
tive nearest neighbors exchange coupling, it is interesting 
to notice that better agreement with experimental data 
can be achieved by proposing a simple dependence of the 
form: 



J(D) = J{oo)e a /°, 



(23) 



in Eq. (ll9l) . This proposal is based on the experimen- 
tal fact that the lattice parameter of metallic nanopar- 
ticles contracts with decreasing particle size in such a 
way that the lattice parameter contraction (Aa/a) is an 
inverse function of the diameter of nanoparticles (25l l33h - 

. Such a lattice contraction is attributed to reduction 
of surface bonds length as a response to surface stress. 
Contraction factors may vary with material and crystal 
orientation. It has been observed experimentally that 
the lattice parameter contracts by 2.4% in 5 nm Ni par- 
ticles [37j |. On the other hand, fee Ni is considered as 
a strong ferromagnet characterized by a less pronounced 
RKKY behavior, exponentially damped, and a faster de- 
cay of the exchange integral with the interatomic dis- 
tance. More concretely, Ni remains ferromagnetic up to 
the 5 nearest neighbors, and within this range of dis- 
tance the exchange integral is essentially a decreasing 
exponential function of interatomic spacing [38, 39]. 

Concerning finite size scaling (FSS) properties, Fig. 5 
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FIG. 5: Log-log plot of the size dependence of the reduced 
temperature for pure nanoparticles having core coordination 
numbers z~S (bcc) (black solid line) and z—12 (fee) (dashed 
line). Gray lines corresponds to the log- log fitting process 
using finite size scaling (FSS) theory (Eq.(f24]|). 



Regarding thermal properties, Fig. 6 shows the tem- 
perature dependence of both the magnetization per site 
and the magnetic susceptibility for pure iron (cc=0, q=0 
and p=l) nanoparticles and for some selected diameters. 
Results for bulk iron are also included for comparison. A 
well-behaved thermal driven ferromagnetic to paramag- 
netic phase transition is observed as well as the shift of 
the critical temperature to low values as the system size 
decreases. The location of the maximum susceptibility 
coincides with that derived from Eq. (116[) . 
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FIG. 6: Temperature dependence of both the magnetization 
per site and magnetic susceptibility for a-Fe nanoparticles 
and bulk iron. 



C. FeMnAl nanoparticles 

In the case of FeMnAl nanoparticles, we have cho- 
sen the stoichiometry Feo.5Mno.1Alo. 4, which, under bulk 
conditions, has been studied by using Mossbauer spec- 
troscopy and magneometric techniques [13] as well as 
from theory [l2| and Monte Carlo simulation [l3j]. As 
is known, the bulk alloy has a Tq close to room temper- 
ature (« 300.fr) in addition to the occurrence of a RSG 
behavior in the low temperature regime. The critical line 
Tc(D) is given by Eq. (|T6"l) . and the corresponding log- log 
plot is shown in Fig. 7. Data have been fitted using finite 
size scaling theory (Eq. (|24j> ). Our best estimate for the 
correlation length critical exponent is v — 0.926 ± 0.004. 
This exponent is still quite similar to that of a 2D Ising 
model, but slightly different from our previous exponent 
for the pure case. This feature is consistent with the Har- 
ris's criterion (46l - |48j for which a different set of critical 
exponents may be expected for diluted and random sys- 
tems having a distribution of exchange integrals. In our 
case, dilution is provided by Al atoms and randomness 
is provided by the random distribution of the atomic el- 
ements in the alloy within the crystalline structure and 
over the entire volume of the nanoparticles. The expo- 
nent is also greater than the computed via Monte Carlo 
simulation (y = 0.79 ± 0.03) of Fe .5Mn .iAlo.4 alloys 
under bulk conditions [l3j. This fact agrees also with 
Fisher's theory for which critical behavior is modified 
when free-edge boundary conditions are considered giv- 
ing rise to a new set of critical exponents different from 
those of the bulk [5J |H|. 

In Fig. 8, we show the temperature dependence of mag- 
netization per bond and magnetic susceptibility as ob- 



tained from Eqs. (fl"4"|) and (fl~5|) respectively. As is ob- 
served, our model predicts a magnetization reduction 
in the low temperature regime, below around 70 K, in 
agreement with zero field cooling measurements for the 
bulk case 



121. 



According to our model, such a reduc- 
tion, which has been attributed to a RSG phenomenol- 
ogy within the ferromagnetic matrix, is still observed for 
nanoparticles. Moreover, the onset of the reentrant phase 
is supported by the low temperature peaks of the mag- 
netic susceptibility. Results reveal also the expected shift 
to lower temperature values of Tc(D) as the particle size 
decreases, in agreement with the critical line obtained 
from Eq.(Tl6|). 



o0.1 



0.01 





Eq.(16) 




Eq.(24) FSS 


' Fe o, Mn o, A1 o, 




.1 , 


, ... 1 



aD (nm) 



10 



30 



FIG. 7: Log-log plot of the size dependence of the reduced 
temperature for Feo.5Mno.1Alo. 4 nanoparticles having core co- 
ordination numbers 2=8 (bec) (circles) using Eq. (|16[) . Solid 
line corresponds to the log-log fitting process for 2=8 using 
finite size scaling (FSS) theory (Eq.((24j)). 

It is well established that Ising spin-glass transitions 
should follow the so-called Almeida-Thouless (AT) line 
[5fjj from which is expected a field dependence of the 
peak temperature obtained from the maximum of 

the magnetization, of the form: 



h oc (1 - T p /T f ) 



3/2 



(25) 



The extrapolation of the AT line at h—0 gives the freez- 
ing SG transition temperature Tf. Agreement of the data 
with the AT line is usually considered as evidence of the 
occurrence of a SG phase, although not concluding [5lj ]. 
Thus, in order to evaluate the properties of the RSG 
phase, we have solved Eg. (1X41 for different low field val- 
ues in Eq. (ri"2"|) from which T p was extracted with an un- 
certainty of ± 1 K. Results are summarized in Fig. 9 
where we plot h 2 / z vs. T p . Two remarkable features are 
observed in this figure. First, our data are in accordance 
with the AT line as one expects for a SG transition at 
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least at low field values. The large plateau observed in the 
magnetization at around T p has been already observed to 
occur experimentally from SQUID and ac susceptibility 
measurements [l4|. Additionally, for large field values, a 
deviation from the A-T line was observed. These results 
are in agreement with those reported by Young et al. [EH 
and references therein where the difficulty of having an 
A-T line for short-range Ising spin glasses at large fields 
was evident. In fact, the existence of a SG ordering in a 
magnetic field is still an open question (53|. Second, the 
freezing temperature is clearly size dependent, i.e. it di- 
minishes as the particle size decreases. This fact implies 
that the SG region in the magnetic phase diagram be- 
comes smaller for nanoparticles exhibiting SG behavior 
within their entire volume and not as a consequence of a 
merely surface effect as it has been proposed in nanopar- 
ticles exhibiting the so-called surface spin-glass-like be- 
havior. These results suggest that SG behavior observed 
in bulk systems is reduced when finite size effects be- 
come important, which could be attributed to a reduc- 
tion in the total number of frustrated spins as the particle 
size becomes smaller. Moreover, in the framework of the 
mean-field approximation, and taking into account that 
the low temperature transition occurs within the ordered 
Fe ferromagnetic matrix, the particle size dependence of 
the freezing temperature can be understood by writing 
Tf = z cS x JMn-Mn + Zefff>JFe-Mn, which means that Tf(D) 
should become proportional to z e g, or at least as 1/D 
according to Ea. tTTTI) . Inset in Fig. 9 reveals that such a 
trend is fulfilled. 

In order to interpret how such a low-temperature 
magnetization reduction takes place, we have performed 
a single-spin flip Metropolis Monte Carlo simulation 
[4ll l54j of Feo.5Mno.1Alo. 4 nanoparticles in the frame- 
work of a nearest-neighbor Ising model. We have used 
free boundary conditions, a maximum of lxlO 5 Monte 
Carlo steps per spin (MCS) and discarded the first 6x 10 4 
MCS for equilibration. Configurational averages over five 
different random atomic realizations were performed. We 
have also employed the same set of competitive parame- 
ters used in the variational approach and numerical val- 
ues of JFe-Fe [H| reproducing the critical temperatures 
under bulk conditions. An example of the simulation re- 
sults for a particle size D—10 is shown in Fig. 10. In 
addition to the total magnetization per site, the corre- 
sponding Fe and Mn contributions are shown separately. 
Both approaches, variational and simulational, predict a 
low temperature magnetization reduction. Monte Carlo 
results allow to conclude that such reduction arises from 
Mn moments for which an increase in the absolute value 
of the corresponding magnetization contribution was ob- 
served. These moments tend to align antiparallel with 
respect to the total magnetization direction ruled by the 
iron matrix according to the negative values of JFe-Mn 
and JMn-Mn- Moreover, such moments are not compen- 
sated and some of them are frustrated. 




Temperature (K) 

FIG. 8: Temperature dependence of magnetization per bond 
and magnetic susceptibility for Feo.5Mno.1Alo. 4 nanoparticles 
and for some selected particle sizes. Data corresponding to 
the bulk case are also included for comparison. 



One of the difficulties of the variational approach is 
that the magnetization obtained in the random diluted 
case is a magnetization per effective bond and not per 
atomic site. This fact gives rise to different values of the 
maximum magnetization (see Fig. 8) in contrast to the 
observed via Monte Carlo, where the maximum value 
of the overall magnetization is close to 0.5 accordingly 
with the Fe atomic concentration, which is practically 
the same independent of the particle size. Finally, re- 
garding the critical behavior along the Tq line, we have 
determined the correlation length critical exponent from 
the maxima of the logarithmic derivative of the magne- 
tization in the vicinity of Tq by assuming the following 
ansatz (45[: 

\ / max 
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FIG. 9: Low field A-T line. Freezing temperature goes down 
as the particle size diminishes. The Tj value for bulk at 
around 62.5 K is relatively close to that reported experimen- 
tally at around 78 K from ac susceptibility measurements for 
Feo.5MnQ.1Alo. 4 bulk alloys [l4| . Inset shows that a T; versus 
1/D linear dependence is followed. 



The log-log plot of the size dependence of the maxi- 
mum values of these derivatives is shown in Fig. 11 
from which our best estimate for the exponent was 
i/=0.71±0.04, very different from that found from the 
variational approach. Our value is however somewhat 
greater than the z^=0.6289±0.0008 value obtained by Fer- 
renberg and Landau (45j | for a 3D Ising model, where 
a high-resolution Monte Carlo study was carried out, 
and somewhat smaller than the z/=0.79±0.03 value ob- 
tained for the same system under bulk conditions f]~3j j . 
These comparisons suggest, contrary to the variational 
approach, that the 3D— s>2D dimensionality crossover does 
not take place, at least for the range of D values we have 
considered in the present study. The difference between 
the correlation length critical exponents, obtained via the 
variational approach and Monte Carlo, is attributed to 
the, so-called, pairs approximation in the variational ap- 
proach where the trial partition function is computed 
by dividing the system only in blocks of one and two 
spins. This ends up in a reduction of the degree of cor- 
relation and consequently in a change of the exponent 
value. Certainly, a more precise estimative of the parti- 
tion function, and consequently of the magnetization and 
Tc, can be achieved by considering the system as formed 
by higher order blocks (four, six, eight spins, etc) [24j |. 
Even though such a calculation is tractable, it turns out 
heavy and very time consuming. This fact constitutes 
the main limitation of the variational model. 
On the other hand, differences respect to the v exponent 
of the pure 3D Ising model can be attributed to the di- 




50 100 150 200 250 300 

Temperature (K) 



FIG. 10: Iron and manganese contributions to the total mag- 
netization per site for Feo.5Mno.1Alo. 4 nanoparticles with di- 
ameter D=10. Inset shows a zoom of the low temperature 
behavior. 



luted character of our system in addition to the disorder 
involved in the distribution of exchange integrals, which 
is consistent with Harris's criterion (40,113], whereas the 
difference with respect to the v exponent of the same sys- 
tem with periodic boundary conditions is consistent with 
Fisher's theory [4(il |. 




10 15 20 25 30 35 



D 



FIG. 11: Log-log plot of the particle diameter dependence 
of the maxima values of the logarithmic derivatives of the 
magnetization to determine the v exponent. 
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IV. CONCLUSIONS 



The critical behavior of ferromagnetic pure and ran- 
dom diluted nanoparticles with competing interactions 
has been addressed. In both cases we have employed the 
free energy variational principle based on the Bogoliubov 
inequality and an Ising model. In the case of random 
diluted nanoparticles, for which we have considered 
the Feo.5Mno.1Alo. 4 system as a case of study, we have 
used, additionally, standard Monte Carlo simulation. In 
order to validate the use of the variational approach in 
nanoparticles, which is carried out for the first time, the 
model was applied to account for the critical behavior 
of pure ferromagnetic nanoparticles on the basis of an 
average nearest neighbors coordination number obtained 
via numerical simulation. Our results allow to conclude 
that the variational model is successful in reproducing 
features of the particle size dependence of the Curie 
temperature for both pure and random diluted parti- 
cles. Comparisons with other theoretical models and 
experimental results for Ni nanostructures were carried 
out in order to evaluate the reliability of the model. A 
better agreement with experimental data was obtained 
if a particle size dependence of the nearest neighbors 
exchange integral is considered consistent with previous 
works where lattice contraction of metallic nanoparticles 
has been observed. 



For random diluted nanoparticles, low temperature 
magnetization reduction was consistent with the same 
type of RSG behavior observed in the bulk counterparts 
in accordance with the Almeida- Thouless line at low 
fields. Such a RSG behavior is attributed to the presence 
of competing interactions, randomness and the aluminum 
dilution effect. A linear dependence of the freezing tem- 
perature with the reciprocal of the particle diameter was 
also obtained indicating that the corresponding region in 
the magnetic phase diagram becomes smaller as the par- 
ticle size diminishes. Concerning critical behavior, data 
obtained by using the variational method were fitted with 
the relationship for critical temperature derived from fi- 
nite size scaling theory and the best estimate for the 
correlation length critical exponent was i/=0.926±0.004. 
Differently from this, a value i/=0.71±0.04 was obtained 
via Monte Carlo. Differences are attributed to the so- 
called pairs-approximation in the variational approach. 
From both approaches, differences in the v exponent of 
Feo.5Mno.1Alo. 4 nanoparticles with respect to that of the 
pure Ising model agree with Harris and Fisher arguments. 
Finally, we want to stress that, even though thermody- 
namical models can be indeed used in the study of nanos- 
tructures [56j j and they can reproduce experimental fea- 
tures, special attention must be paid regarding critical 
behavior depending on the approximations of the model. 
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